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Abstract. The use of polarimetric techniques is nowadays widespread among 
solar and stellar astronomers. However, notwithstanding the recommandations 
that have often been made about the publication of polarimetric results in the 
astronomical literature, we are still far from having a standard protocol on which 
to conform. In this paper we review the basic definitions and the physical sig- 
nificance of the Stokes parameters, and we propose a standardization of the 
measurement of polarized radiation. 



1. Introduction 

The polarization properties of a radiation beam can be described in several 
different ways but have ultimately to rely on the specification of four independent 
quantities. Particularly suitable from this point of view are the so-called Stokes 
parame ters w hich were introduced in the scientific literature by George Gabriel 
Stokes and which are commonly referred to, in modern notations, with 

the symbols /, Q, U, and V. It is possible to give two different definitions for the 
Stokes parameters of a beam of electromagnetic radiation propagating along a 
given direction. According to the first one, the Stokes parameters are expressed 
as statistical averages of bilinear products of the components of the electric field 
along two perpendicular axes, both perpendicular to the direction of propagation 
of the radiation (see Sect. 3). This basic definition is also used in practice by 
radio- astronomers for expressing their measurements of the Stokes parameters. 
The second definition, mostly used by optical astronomers, is an operational one 
that makes use of the concept of ideal filters (see Sect. 7). Obviousy, the two 
definitions have to be consistent. 

2. The choice of the reference system 

Before introducing these two definitions, it is fundamental to point out that 
both of them need the preliminary choice of a reference direction pertaining to 
the plane perpendicular to the direction of propagation of the radiation. To 
introduce the first definition of the Stokes parameters (the one based on bilin- 
ear products of the components of the electric field), we refer to a right-handed 
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coordinate system, (x, y, z), with the z-axis directed along the direction of prop- 
agation and the x-axis directed along the reference direction. In astronomical 
observations, the reference direction can be thought of as the great circle per- 
taining to the celestial sphere and passing through the observed object. Though 
the choice of the great circle is, in principle, arbitrary, it has become custom- 
ary, at least in night-time astronomy, to choose the reference direction as the 
celestial meridian passing through the observed object. In solar observations 
this convention is not always followed since there are, in general, more physical 
directions to choose, like the tangent to the solar limb for prominence obser- 
vations, or the direction passing through the observed point and the center of 
the solar disk for sunspots observations. In polarimetric measurements of other 
solar system objects (e.g., asteroids) it is common to adopt as reference direction 
the great circle passing through the object itself and the Sun, or, better choice, 
the perpendicular to it (see Sect. 6). 



3. The basic definition 

According to Maxwell equations, the electric field vector of the radiation beam 
lies in the {x, y) plane so that it is described by only the two components Ex and 
Ey. In a fixed point of space along the direction of propagation, at the entrance of 
the telescope for instance, these two components of the electric field are described 
by the functions of time, Ex{t) and Ey(t), for which we can introduce the Fourier 
transforms, Sxiijj) and Sy{u!), according to the usual definitions 

£A^) = /_°°^i?.(t)e--*di 

The Stokes parameters of the radiation beam at the angular frequency lu are 
defined by the expressions 

liu) = k[{£x{ojr£x{oj)) + {£y{u;y£y{u;))] 

Qico) = k[{£x{ojr£x{u^))-{£y{u;)*£y{u;))] 

Uiio) = k[{£xiu;)*£yiu;)) + {£yiu;)*£xiio))] 

V{u;) = ik[{£x{ujr£y{u;))-{£y{u;r£x{iv))] 

where fc is a positive constant, the symbol * means complex conjugate, and 
the symbol (• • •) means a statistical average which is implicit in the process of 
measurement. 

There are some risks of ambiguity in the definition of Stokes parameters. 
One of course is the choice of the handedness of the reference system. We have 
adopted a right handed reference system, but adopting a left handed system will 
change the signs of U and V. 

There is a much more subtle ambiguity arising from the definition of Fourier 
transform. In our definitions of Eq. we have implicitely adopted for the 
Fourier transform the representation of Eq. Introducing the alternative 

definition with the substitution 
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namely, defining 

one obviously has, being Ex{t) and Ey{t) real, 

so that the sign of the V Stokes parameter changes when the definition of the 
Fourier transform is changed from Eq. (1) to Eq. (3). In practice, this means 
that when adopting Eq. (2) for the definition of the Stokes parameters, one has 
to be aware that such definition implies a convention about the choice of positive 
or negative time in the definition of the Fourier transform. 

Finally, we want to comment on the use of an arbitrary constant in front 
of the bilinear expressions that defines the Stokes parameters. As a matter of 
fact, we always use relative quantities, that is, we normalise Q, U, and V to 
the intensity, and we consider the reduced Stokes parameters Pq, Pu, and Py 
defined by 

Pq = Q/l 

Pu = U/I (4) 
Pv = V/I 

so the value of the constant k can be left undefined. 



4. The physical significance of the Stokes parameters 

According to the definition given by Eqs. @, Stokes Q is the difference between 
the amount of photons whose electric field oscillates along the reference direc- 
tion and the amount of photons whose electric field oscillates in the direction 
perpendicular to it; Stokes U is the difference between the amount of photons 
whose electric field oscillates at 45° and the amount of photons whose electric 
field oscillates at 135° with respect to the reference direction (the angles are 
reckoned counterclockwise from the reference direction looking at the source - 
see Fig. 1). 

Stokes V is given by the so called right handed circular polarization minus 
the left handed circular polarization, that are so defined: in a fixed point of 
space, the tip of the electric field vector carried by a beam having positive 
circular polarization rotates clockwise, as seen by an observer looking at the 
source of radiation. Viceversa, the tip of the electric field vector carried by a 
beam having negative circular polarization rotates counterclockwise, as seen by 
an observer looking at the source (see Fig. 

It is possible to show that Stokes parameters give information about the 
ellipse drawn by the tip of the electric vector in the plane perpendicular to the di- 
rection of wave propagation. For the sake of simplicity, one can consider the case 
of a monochromatic wave, which is completely polarized. More generally, we re- 
mind that any radiation beam can be considered as the incoherent superposition 
of an unpolarized beam and a beam which is totally polarized (Chandrasekhar 
theorem). The tip of the electric field associated with this beam draws, in a 
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Figure 1. Geometrical visualization of the definition of Stokes parameters. 
The observer is looking toward the source. 



fixed point of space, an ellipse. The major and minor semi-axes of this ellipse, a 
and b, (expressed in non dimensional units), and the tilt angle, 6, that the major 
semi-axis forms with the reference direction x (reckoned counterclockwise -for 
an observer looking at the source- from the reference direction) are connected 
to the reduced Stokes parameters of the polarized component by the equations 



a 
b 

Pu/Pq 



1 - Pv.)l/2 + (1 
l|(l_p^)l/2_(l_ 

tan (2(9) 



hPy)l/2 

Pv)'/'\ 



(5) 



where \z\ means absolute value of z. Stokes V is connected to the eccentricity of 
the polarization ellipse, and Stokes Q and U tell us how the polarization ellipse 
is oriented. If linear polarization is zero, the polarization ellipse degenerates in 
a circle. If circular polarization is zero, the polarization ellipse degenerates in a 
segment. 



5. Linear polarization and position angle 

Linear polarization is often expressed also in terms of Pl and 6, where Pl is the 
fraction of linearly polarized radiation, and 9 is the angle of maximum polariza- 
tion, i.e., the angle that the major axis of the polarization ellipse forms with the 
X-axis of the reference system reckoned with the usual convention: 

^ = ^/W^ 

Pq = Plcos(20) (6) 
Pu = PLSin(2^) 

The use of Pl and 9 is widespread in astronomy because it allows the graphical 
superposition of linear polarization results to the image of the observed object. 
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Note that from Eq. Q one may be tempted to calculate the position angle 9 
as (1/2) arctan (P[//Pq). This is incorrect, as can be seen with the following 
example: consider the case oi 9 = 67.5° i.e., a polarization ellipse with its 
major axis tilted at an angle of 67.5° with respect to the reference x-axis. From 
Eq. we obtain Pq = —Pi^j^fl and P\j = Pi/\/2. By blindly computating 
(1/2) arctan (P(7/Pq) one will find an ellipse tilted at —22.5°, i.e., perpendicular 
to the original, correct one. The reason for that inconsistency is that arctan is 
a function defined between —90° and 90°. The correct expression for obtaining 
the position angle (apart from inessential multiples of vr) is the following one: 

9 = +^7r if = and Pu > 

9 = -Itt if i^Q = and Pu < 

9 = \ arctan (^) if Pq > ( 

9 = \ arctan j + f if Pq < 

Alternatively, one can substitute Eqs. (O by the unique formula, suggested by 
V. Bommier (private communication) 

9 = ^ sign(Pt/) arccos ( ^J'^ | . (8) 




6. Transforming Stokes parameters into a new reference system 

A further point concerns the transformation law for the Stokes parameters when 
changing the reference direction. If the new reference direction is obtained from 
the old one by a counterclockwise (looking at the source) rotation of an angle Xi 
by means of Eq. @ one can easily show that the new reduced Stokes parameters, 
{Pq, PIj, Py)j are connected to the old ones, {Pq,Pi/,Pv), by the equations 

P^ = cos(2x)PQ+sin(2x)Pc/ 

PIj = -sin(2x)PQ + cos(2x)Pc/ (9) 
P(^ = Pv 

Equations © show that the linear polarization Stokes parameters, Q and U, 
transform between each other when the reference direction is changed, whence 
the importance, as already remarked, of clearly specifying the chosen reference 
direction when publishing a result concerning linear polarization. The same 
equations also show the invariance of the quantity Pl under a rotation of the 
reference direction, and the obvious fact that a rotation of the reference direction 
of 180° leaves all the Stokes parameters invariant. This last fact means that when 
choosing the reference direction it is not necessary to specify "an arrow" on the 
direction itself. If, for instance, one chooses as reference direction the celestial 
meridian passing through the observed object, it does not matter whether the 
X-axis of the system (x, y, z), on which the electric field components are defined, 
is directed towards the North or the South pole. 

It should be noted that Eq. @ says that the position angle 9 represents the 
angle by which one should rotate the reference direction (in the counterclockwise 
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direction, looking at the source), in order to end up with a positive Pq value 
(which indeed turns up to be equal to Pl) and a null Pu value. 

As an example of the use of Eq. ©, let us consider the polarization of 
a solar system object measured in a system that has the reference direction 
(x-axis) aligned to the celestial meridian passing through the observed object. 
We want to transform Stokes parameters into a system that has the reference 
direction perpendicular to the great circle passing through the object and the 
Sun. We have to calculate Pq, PIj, Py via Eqs. ©, setting x = ^ + 7r/2, where 
$ is the angle between between the celestial meridian passing through the object 
and the great circle passing through the Sun and the object. The angle $ is 
related to the coordinates of the object and of the Sun through the relationship 

sin5TCOs(a0 — ax) = tan^© cos^t — sin(a0 — ot) (10) 

tan 9 

where (q©,*^©), and (otj^^t) are the right ascension and declination of the Sun, 
and of the observed object, respectively. Alternatively, one can directly calculate 



p/ _ B'^-C^ p 2BC p 

- ~ ^^^^ 



p/ _ 2BC p I B^-C^ p 



where 



A = sin 5q sin 5t + cos 5q cos 5t cos (a© — ax) 

B = cos ^0 sin(a0 — ax) (12) 

C = sin ^0 cos Jt — cos 5q sin 5^ cos(a0 — ax) • 

This transformation is meaningfull in polarimetry of solar system objects be- 
cause in the new reference system. Stokes Q represents the flux perpendicular 
to the plane Sun-Object-Earth (commonly referred to as the scattering plane) 
minus the flux parallel to that plane. 



7. Operational definition 

In optical polarimetry, it has become customary to give an alternative definition 
of the Stokes parameters in terms of ideal filters for linear and circular polar- 
ization (Shurcliff 1962). This requires the preliminary definitions of ideal filters 
for linear and circular polarization. 

7.1. An ideal filter for linear polarization 

An ideal filter for linear polarization (called a polarizer) is a device that can be 
interposed along a radiation beam and which, by definition, is totally transparent 
to the component of the electric field along a given direction, perpendicular to 
the direction of propagation (the transmission axis of the polarizer) , and totally 
opaque to the component of the electric field in the orthogonal direction. 

7.2. An ideal filter for circular polarization 

An ideal filter allowing for positive circular polarization is a device that is totally 
transparent to a radiation beam whose electric field component along the y-axis 
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has, in a fixed point of space, a phase lag of 90° with respect to the component 
along the x-axis, and is totally opaque to a radiation beam which presents a 
phase lag of —90° (or a phase anticipation of 90°) between the same components. 
An ideal filter allowing for negative circular polarization acts in the opposite 
way. Note that according to the previous statement, and consistently with our 
choice of the right-handed coordinate system {x,y,z), right and left circular 
polarization are those defined in Sect. 4. 

An ideal filter allowing for positive circular polarization can be realized by 
the combination of an ideal quarter-wave plate followed by an ideal polarizer 
whose transmission axis is rotated counterclockwise (looking at the source of 
radiation) by an angle of 45° with respect to the direction of the fast-axis of 
the plate. ^ An ideal filter allowing for negative circular polarization can be 
realized by the combination of an ideal quarter-wave plate followed by an ideal 
polarizer whose transmission axis is rotated clockwise (looking at the source of 
radiation) by an angle of 45° with respect to the direction of the fast-axis of 
the plate. Note that from the practical point of view it is common to consider 
a retarder waveplate free to rotate, followed by a fixed linear polarizer, the 
transmission axis of which is aligned to the reference direction x. Obviously, 
in this configuration, the ideal filter allowing for positive circular polarization 
is obtained setting the fast-axis of the retarder waveplate at an angle of —45° 
(i.e., 45° clockwise) with respect to the reference direction, looking at the source 
of radiation, and the ideal filter allowing for negative circular polarization is 
obtained setting the fast-axis of the retarder waveplate at an angle of 45° with 
respect to the reference direction. 



7.3. The operational definition of Stokes parameters 

The operational definition of the Stokes parameters of a radiation beam can be 
given in terms of the ideal filters introduced above and of an ideal detector, 
capable of producing a signal proportional to the electromagnetic energy falling 
on its acceptance area in a given interval of time and in a given interval of 
angular frequency (v, u + di^). We denote by S the signal of the detector with no 
filter interposed and by 5*0, 545, 5*90, and 5*135 the signals of the detector with 
an ideal filter for linear polarization interposed with the transmission axis set, 
respectively, at 0°, 45°, 90°, and 135° with respect to the reference direction, all 
the angles being reckoned counterclockwise looking at the source. Finally, we 
denote by S-^- and S- the signals of the detector with an ideal filter interposed, 
allowing, respectively, for positive and negative circular polarization. The Stokes 
parameters are defined by the equations 



^Wavoplatos arc optical elements with two principal axes, one called the fast-axis, and the other 
one the slow-axis, characterised by two difltercnt refractive indices. The refractive index of 
the fast-axis is smaller than the refractive index of the slow-axis, so that linearly polarized 
light with its electric vector parallel to the fast-axis travels faster than linearly polarized light 
with its electric vector parallel to the slow-axis. A quarter-wave plate produces a n/2 phase 
retardation between the components of the electric field along the fast and slow axes. 
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where k' is a positive constant that turns out to be inessential when considering 
the ratios Pq, Pu, and Py of Eq. (jlj). 

8. Conclusions 

It can be easily shown that the two definitions of the Stokes parameters (the one 
based on the electric field components of Eq. ^ in Sect. 3, and the operational 
one, based on the concept of ideal filters of Eq. ((T^ in Sect. 7.3) are completely 
equivalent in as far as the fractional Stokes parameters are concerned. The op- 
erational definition is the same as the one proposed by Shurcliff (1962) and both 
definitions agree with those of Landi Degl'Innocenti and Landolfi (2004). We 
believe that the astronomical community (including radio-astronomers) would 
certainly profit in conforming to this standard. Finally, although each reference 
direction is as good as any other in as far as its choice is clearly specified, we 
propose that, except for observations of objects belonging to the Solar System, 
the astronomical community should conform to the standard of adopting as ref- 
erence direction the celestial meridian passing through the observed object. For 
solar system objects other than the Sun, we propose to adopt as reference di- 
rection the perpendicular to the great circle passing through the object and the 
Sun. 
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